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ABSTRACT 

Beginning with a simple model for the growth of structure, I consider the dissipation- 
less evolution of a MOND-dominatcd region in an expanding Universe by means of 
a spherically symmetric N-body code. I demonstrate that the final virialized objects 
resemble elliptical galaxies with well-defined relationships between the mass, radius, 
and velocity dispersion. These calculations suggest that, in the context of MOND, 
massive elliptical galaxies may be formed early (z > 10) as a result of monolithic 
dissipationless collapse. Then I reconsider the classic argument that a galaxy of stars 
results from cooling and fragmentation of a gas cloud on a time scale shorter than 
that of dynamical collapse. Qualitatively, the results are similar to that of the tradi- 
tional picture; moreover, the existence, in MOND, of a density-temperature relation 
for virialized, near isothermal objects as well as a mass-temperature relation implies 
that there is a definite limit to the mass of a gas cloud where this condition can be 
met- an upper limit corresponding to that of presently observed massive galaxies. 



1 INTRODUCTION 



Here I reconsider the problem of galaxy formation, dissipa- 
tionless and dissipational, in the context of modified New- 
tonian dynamics, or MOND (Milgrom 1983). With New- 
tonian gravity, dissipational processes do not appear to be 
necessary to explain the basic structural properties of el- 
liptical galaxies. For example, van Albada (1982) demon- 
strated that the observed surface brightness profile of ellip- 
tical galaxies, the r 1 ^ 4 law (de Vaucoleurs 1948), develops 
naturally in Newtonian N-body simulations of the collapse 
of bound systems with initially inhomogeneous density dis- 
tributions. On the other hand, dissipation does seem to be 
necessary to explain the typical densities of baryonic matter 
in both disk and elliptical galaxies (Binney 1977) ; pure grav- 
itational collapse in an expanding Universe produces insuffi- 
cient concentration of visible matter. Moreover, cooling and 
fragmentation are certainly necessary for the formation of 
stars which are the principal component of galaxies. Thus, 
the formation of galaxies could, in some sense, be viewed 
as dissipationless if cooling of a primordial gas cloud and 
fragmentation down to the level of stellar mass objects oc- 
curs on a time scale short compared to that of dynamical 
collapse. This is a point first made by Hoyle (1953) to ex- 
plain the mass scale of galaxies, and developed further by a 
number of authors (Silk 1977, Binney 1977, Rees & Ostriker 

1977, White & Rees 1978). The argument involving these 
two competing timescales has become an essential ingredi- 
ent of the standard cosmology in which the matter budget 
of the Universe is dominated by a hypothetical pressureless, 
collisionless fluid, cold dark matter, or CDM (White & Rees 

1978, Blumenthal et al. 1984). 

In the context of the MOND paradigm, galaxy for- 
mation must occur without the assistance of CDM. More- 
over, recent numerical work has indicated that dissipation- 



less merging of galaxies proceeds much more slowly with 
MOND than when luminous galaxies are assumed to be sur- 
rounded by massive dark halos (Nipoti, Londrillo & Ciotti 
2007b); therefore, mergers would probably play a less impor- 
tant role in galaxy formation, particularly the formation of 
ellipticals, than in the current conventional picture. But the 
primary roadblock to the consideration of galaxy formation 
with MOND is that there is not yet a generally accepted or 
standard cosmological context for structure formation, even 
though relativistic extensions have been described in the lit- 
erature (e.g., Bekenstein 2004, Sanders 2005, Zlosnik, Fer- 
reira & Starkman 2006). Therefore, I begin with the simple 
assumptions that MOND only applies to peculiar acceler- 
ations in an expanding Universe and that the acceleration 
constant, ao, does not evolve with cosmic time. 

Such assumptions would appear to be consistent with 
the current relativistic versions of MOND, where the non- 
Newtonian force is mediated by a long-range scalar field with 
non-standard Lagrangian: there is no MOND in the absence 
of density variations. Moreover, these assumptions underly 
the heuristic models of structure formation in a MONDian 
universe considered previously by myself (Sanders 2001) and 
by Nusser (2002). In particular, Nusser noted that small fluc- 
tuations grow rapidly (oc (z + 1)~ 2 ), and that this leads to a 
unacceptably clumpy Universe at the present epoch- essen- 
tially independent of the initial amplitude of fluctuations. 
I postpone consideration of this problem to a later discus- 
sion on taming the growth of large-scale structure formation 
and assume here that galaxy scale fluctuations in a baryonic 
Universe grow by the application of the MOND formula to 
peculiar accelerations. If so, then it is easily demonstrated 
that galaxy mass objects re-collapse early (10 < z < 25) and 
that spherically symmetric dissipationless collapse leads to 
final virialized objects having roughly the properties of ob- 
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served elliptical galaxies. In particular, both the magnitude 
and form of the observed surface density distribution -the de 
Vaucoleurs law- and the velocity dispersion-baryonic mass 
relation- the Faber- Jackson law (Faber & Jackson 1976)- 
are reproduced for bound objects with masses ranging up to 
10 13 M Q . 

All of this, however, begs the question of why galaxies, 
as self-gravitating ensembles of stars, have an upper mass 
limit to their baryonic content of about 10 12 Mq. More mas- 
sive self-gravitating objects exist, but they are bound groups 
or clusters of galaxies rather than single entities consisting 
of stars. To address this question we must again consider 
the processes of radiative cooling and fragmentation on a 
time scale short compared to that of gravitational collapse. 
In this respect, there is a great advantage afforded by mod- 
ified dynamics: the presence of an additional physical con- 
stant with units of acceleration (ao ~ 1CF 8 cm/s 2 ) in the 
structure equation strongly constrains the properties of any 
self-gravitating object as we see in this and previous dissi- 
pationless collapse calculations (Nipoti, Londrillo & Ciotti 
2007a). Such an object with a velocity dispersion of a few 
hundred km/s will have a mass of 10 11 M as well as a 
characteristic size and density. In other words, there is less 
arbitrariness in specifying initial properties of a pre-galactic 
cloud. 

In the context of MOND the virial mass, the dynamical 
timescale, and the characteristic density of a self-gravitating 
isothermal cloud depend only upon velocity dispersion. This 
redefines the curve, in the temperature-density plane, de- 
scribing the condition that dynamical and cooling timescales 
are equal. As in the standard model, this curve neatly 
separates those objects where dissipation has led to the 
formation of stars (galaxies) from those where it has not 
(groups and clusters). Moreover, the presence of a density- 
temperature relation for self-gravitating isothermal clouds 
leads to definite upper limit to the mass of objects in which 
the cooling time is less than the dynamical time- again a 
characteristic mass of about 10 12 M®. Therefore, we see that 
many of the successes of the standard model for galaxy for- 
mation carry over to MOND- in particular the existence 
of an upper limit to systems consisting of stars. Moreover, 
the location of near isothermal systems in the temperature- 
density plane- systems ranging from globular clusters to 
clusters of galaxies- plane becomes understandable in terms 
of MOND. 

In the following section I review the heuristic model for 
the initial growth of galaxy scale fluctuations in an expand- 
ing Universe. In section 3 I numerically follow the dissipa- 
tionless evolution of galaxy-scale spherical regions which col- 
lapse out of the Hubble flow. Even though dissipationless, 
these calculations demonstrate the scaling relations which 
primordial gas clouds may obey. In section 4 1 apply these re- 
lations to reconsider the old argument on collapse and frag- 
mentation in the context of MOND, and in the final section 
I compare MOND with the standard model and summarise. 



2 THE GROWTH OF SMALL FLUCTUATIONS 
IN A MONDIAN UNIVERSE 

MOND in its original form (Milgrom 1983) is described by 
the relation between the true gravitational acceleration, g, 



and the Newtonian gravitational acceleration, g^— 

gn(g/a )=g N (1) 

where ao is the fundamental acceleration parameter (~ 10~ 8 
cm/s 2 ) and /j, is the function which interpolates between the 
Newtonian regime (g > ao, = 1) and the MOND regime 
(g < ao, fi = g/ao). MOND, as a modification of Newtonian 
gravity, may also be expressed by the modified field equation 
of Bekenstein & Milgrom (1984): 

V ■ [fj,{V4>/ao)V4>] = 4-irGp. (2) 

This formulation is conservative but difficult to solve for ar- 
bitrary mass distributions; in the case of spherical symmetry 
it reduces to the simple expression above (eq. 1) where gjv 
is determined by the usual Poisson equation. 

Early efforts to describe a MONDian cosmology in the 
absence of a relativistic theory applied the simple MOND 
equation (eq. 1) to an expanding spherical region analo- 
gously to the Newtonian derivation of the Friedmann equa- 
tion (Felten 1984). There it was immediately realized that 
the evolution depended upon the physical size of the region 
(it was not possible to define a dimensionless scale factor) 
and that any expanding region would eventually re-collapse 
regardless of its initial density and expansion velocity. How- 
ever, at early epochs, the region within which the Hubble 
deceleration is less than ao, is much smaller than the horizon 
scale; therefore, if ao is independent of cosmic time, it may 
be possible that the Universe as a whole is Friedmannian 
while increasingly larger sub-regions become MONDian and 
re-collapse (Sanders 1998). Although this would lead natu- 
rally to a scenario of hierarchal structure formation, it is not 
possible to identify the centres of re-collapse; i.e., primordial 
density fluctuations play no role. 

A more plausible scenario would be one in which the 
MOND formula is applied to the peculiar accelerations de- 
veloping from density fluctuations rather than to the Hub- 
ble expansion as a whole. This was the central idea be- 
hind heuristic models for structure formation (Sanders 2001, 
Nusser 2002) and for the early collapse of low-mass gas 
clouds (Stachniewicz & Kutshcra 2005). This assumption is 
consistent with, but not a necessary consequence of, cur- 
rent relativistic extensions of MOND (for example, TeVeS, 
Bekenstein 2004). Its validity depends upon the form of the 
free function of the theory (effectively the MOND interpo- 
lating function) in the cosmological regime where the cosmic 
time derivative of the scalar field dominates the scalar field 
invariant. But since the free function does not follow from 
any more fundamental considerations, assumptions about 
its form, at present, are no less ad hoc than applying the 
MOND prescription directly to peculiar accelerations. 

With the additional ansatz that the MOND acceleration 
parameter ao is does not vary with cosmological time (also 
not necessarily true in relativistic extensions such as the 
biscalar variant of TeVeS, Sanders 2005), the equation for 
the growth of small fluctuations ((S = 8p/p) becomes: 

S + 2 *6+-8 = ^ (3) 

X x xX c 

Here x is the dimensionless scale factor in terms of the 
present scale factor (x(to) = 1), time is in units of the Hub- 
ble time (1/Ho), A c is the co-moving scale of the perturba- 
tion, and gi is the peculiar gravitational acceleration related 
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to the Newtonian peculiar acceleration 

9p = ^A c 5 (4) 

by the MOND formula (eq. 1). Here I have assumed that 
fluctuations grow primarily during the period when non- 
relativistic matter with density parameter O m dominates 
the mass-energy budget of the Universe. Below 1 will as- 
sume, consistent with MOND, that this matter is essentially 
baryonic. 

In the high acceleration limit (gi > do) we recover the 
usual linear Newtonian expression for the growth of small 
fluctuations in an expanding medium (gi — g p ). However, 
when gi < ao, then the MOND limit applies: 



3a; 2 



(5) 



where m is the Hubble radius, c/Ho, and the MOND accel- 
eration parameter is written as ao = f m cHo with f m ~ 1/7. 
In this case, the equation for the growth of small fluctua- 
tions becomes non-linear (g± oc S 1 ^ 2 ) and dependent upon 
the co-moving or mass scale. 

In the MOND limit a simple power law solution to this 
equation exists in the matter dominated regime; i.e., 8 = 
At a with q = 4/3. In other words, unlike the Newtonian 
case where 8 oc x, in the MOND regime 8 oc x 2 (Nusser 
2002). In terms of the scale factor, the solution is 



8 - —f — — 

m 



(6) 



In fact, as may occur in non-linear problems, this solution is 
an attractor on a space of solutions; MOND drives structure 
growth toward x 2 quite independently of the initial condi- 
tions. 

This is shown in Fig. 1, where 8 is plotted as a function 
of scale factor. The cosmological background here is of the 
usual Friedmann form for a Universe without dark matter 
(f2 m = 0.04) but with the standard radiation content and 
zero spatial curvature, i.e., Q\ = 0.96. Note that here I do 
not consider the possible presence of 2 eV neutrinos; this 
would change details- on the scale of galaxies, adding an 
additional contribution to Hubble expansion- but not the 
essence of the argument. In eq. 1 I have taken the form 
of n to be that advocated by Zhao & Famaey (2006), i.e., 
fj,(x) = x/(l + x). The heavy solid curve is the pure MOND 
power law solution (eq. 6) for a co-moving scale of 1.56 Mpc 
corresponding to a mass of 10 11 Mq, and the lighter solid 
curves are numerical solutions of eq. 3 with various initial 
values for <5 at the epoch of matter decoupling (z « 1400). 
We see that the numerical solutions rapidly converge to the 
power law solution. 

It is also evident that the fluctuations grow to unity on 
a at relatively early epoch: 



x x = 2.6 



(7) 



- fm r H . 

corresponding to a redshift of zi = 
13.8(H) 11 M©/M)5fi m _ 3 « 35 for the galaxy scale 
fluctuation. In this scenario we would expect massive 
galaxies to form early (z > 10). 

Given that the power spectrum of fluctuations is re- 
lated to the amplitude as P(k)k 3 oc 8 2 this gives a final 



t>i) 

c 




Figure 1. The amplitude of a galaxy scale density fluctuation 
(A c = 1.56 Mpc) as a function of scale factor on a log-log plot. 
The heavy solid line is eq. 6, the power-law solution of eq. 3, and 
the light solid curves are the numerical solutions of the growth 
equation (eq. 3) in the low density (S7 m = 0.04), vacuum energy 
dominated cosmology for initial values at decoupling: <5o = 1.8 X 
10~ 7 , 1.8 XlO -5 , 1.8 xlO -3 . The dashed curve is the evolution of 
over-density of the expanding sphere determined by the numerical 
N-body program; i.e., when S > 1 cq. 6 is no longer valid. 



power spectrum of P(k) oc fc _1 as noted by Nusser (2002). 
However, also as noted by Nusser, eq. 6 would certainly im- 
ply that the amplitude of fluctuations on large scale is far 
too large to be consistent with observations. Cosmology can 
intervene to tame the growth of larger scale structure; essen- 
tially, vacuum energy (or curvature) moderates the growth 
when it begins to dominate the cosmic expansion, but this 
will be considered in a later paper. For now, we assume 
that eq. 6 applies to galaxy scale fluctuations. Of course, the 
equation is only valid when 8 < 1; thus below we consider the 
dissipationless evolution of collapsing spheres with MOND; 
i.e., only (modified) gravity affects the evolution. The un- 
derlying implicit assumption is that the baryonic fluid, pre- 
sumably a gas cloud at the beginning of re-collapse, rapidly 
cools, fragments and forms stars during the initial collapse. 
The condition for this to occur will be considered in Section 



3 THE EXPANSION AND RE-COLLAPSE OF 
MOND DOMINATED FLUCTUATIONS 

Fully three-dimensional dissipationless collapse calculations 
have recently been carried out by Nipoti et al. (2007a). 
Starting with an inflated Plummer sphere (zero initial ki- 
netic energy) they follow the collapse and formation of viri- 
alized objects with varying central concentration, ranging 
from deep MOND to Newtonian. Here I use a spherically 
symmetric N-body code of the form originally developed by 
Henon (1964) in order to follow the evolution of an initially 
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expanding over-dense spherical region beyond 5=1. Al- 
though the calculations of Nipoti et al. would provide more 
realistic models of actual galaxies, the goal here is to con- 
sider the general properties of objects which might actually 
condense out of the Hubble flow. 

The radial motion of spherical shell i at radius r» is 
determined by numerically solving the equation of motion: 



d 2 n 

dt 2 



= -9* + 



(8) 



where the gravitational force at the i shell, gi, is given by 
two parts: g t = g H + g 1 where g H is the usual Hubble decel- 
eration and gi, as above, is the peculiar deceleration result- 
ing from the over-density and given by the MOND formula 
(eq.. 1), again with the Zhao-Famaey interpolating function. 
Whenever the peculiar acceleration exceeds the Hubble ac- 
celeration by a factor of two, I assume that the region is 
decoupled from the Hubble flow and then follow its evolu- 
tion as an isolated object. The final term in eq. 8 is the 
centrifugal acceleration with ji being the specific angular 
momentum of the i th shell; i.e., ji = v t ri where v t is the 
tangential velocity of particles comprising that shell. Note 
that there is no systematic rotation; vt of particles in a shell 
is assumed to be distributed uniformly in all tangential di- 
rections. Again, because of the perfect spherical symmetry 
the solution of eq. 1 is equivalent to solving the Bekenstein- 
Milgrom modified Poisson equation. 

The sphere is initially homogeneous and partaking in 
the uniform Hubble expansion perturbed by its peculiar 
velocity. An initial over-density of 5p/p — 0.03 is simu- 
lated by taking the sphere to be initially slightly smaller 
(Sr/r = 0.01) than it would be in a uniform Universe. The 
initial conditions for the expanding sphere then are taken 
from the numerical solution for 5 described above (as in Fig. 
1) in the case where the fluctuation amplitude at decoupling 
is So — 1.8 x 10~ 7 . For the sphere with mass 10 11 M©, when 
5 = 0.03 we find U = 0.00278 in units of the Hubble time, 
corresponding to Xi = 0.00968 which yields an initial radius 
of 15.1 kpc and an expansion velocity of 254 km/s (this is 
all for the adopted cosmological model). The expansion is 
taken to be linear (= rH) out to the maximum radius. The 
initial tangential velocity on all shells is assumed to be 20- 
30 km/s in order to prevent a collapse which is too violent 
(the results are qualitatively insensitive to the actual value 
of vt)- With MOND, the over-dense sphere will inevitably 
re-collapse after reaching a maximum radius about 3 times 
larger than the initial radius. The re-collapse occurs after a 
time interval comparable to the age of the Universe at the 
epoch when the sphere decouples from the Hubble flow. Of 
course, the assumptions of uniformity and spherical symme- 
try are enormous simplifications; since smaller regions re- 
collapse earlier the region forming a galaxy would consist of 
numerous collapsing or virialized sub-components. 

This idealised evolution is illustrated in Fig. 2 for the 
10 11 Mq sphere simulated by 800 spherical shells. This shows 
the logarithm of the radius of four different shells as a func- 
tion of time; the largest shell is the outermost shell. It is ev- 
ident that the shells oscillate with different frequencies and 
inter-penetrate. Due to the phase mixing, the sphere will 
eventually come into an equilibrium state. The over-density 
is calculated by comparing the average density within in the 
outermost shell to that of the homogeneous Universe at that 
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Figure 2. The log radius of four characteristic shells in a 10 11 Mq 
sphere, as a function of time. Initially the entire sphere is uni- 
formly expanding and, in a Newtonian context, would be un- 
bound; i.e., the sphere would never re-collapse. But with modified 
dynamics, all shells will eventually re-collapse. After the initial 
re-collapse, the shells oscillate at different frequencies, phase mix, 
and come into equilibrium. 



epoch. This is shown by the dashed curve in Fig. 1 where we 
see that evolution of 5 follows that of the numerical solution 
of eq. 3 until 5 w 1 and then diverges. 

The approach to final equilibrium can be seen by consid- 
ering the virial ratio (2T/V) as a function of time. In modi- 
fied dynamics the gravitational potential is, properly speak- 
ing, not defined; an infinite energy is required to move any 
shell to infinity. None-the-less it is possible to define a virial 
relation which, in equilibrium, takes the form 2T — V — 
where T is the total kinetic energy as usual and, in spherical 
symmetry, 



7 



V = — 47r / p(r)g(r)r dr 



(9) 



(Romatka 1992, Gerhard & Spergel 1994). Here, g(r) is the 
modified gravitational acceleration given by eq. 1. For the 
system of spherical shells this becomes V = rirriig(ri) 
over all shells. As in Newtonian dynamics the condition 
2T/V = 1 corresponds to equilibrium, but, unlike Newto- 
nian dynamics, if 2T/V > 2 the system remains gravita- 
tionally bound and will re-collapse. 

The virial ratio as a function of time for the 10 11 M© 
sphere is shown in Fig. 3. Initially, there is too much kinetic 
energy in expansion as would be the case in a low £7 universe; 
the sphere would never re-collapse in the context of Newto- 
nian dynamics. But with modified dynamics re-collapse and 
mixing do occur and after roughly three or four dynamical 
time-scales (~ 10 9 years), the virial ratio approaches one. 
Fairly significant (~ 10%) oscillations continue for another 
three to four dynamical timescale- somewhat longer than 
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Figure 3. The virial ratio (2T/V) as a function of time for the 
10 11 Mo spherical protogalaxy of Fig. 7 (2T/V = 1 in equilib- 
rium). This illustrates the rapid approach to virial equilibrium 
after entering the MOND regime and re-collapsing. The spherical 
galaxy is in place as a virialized system 2 X 10 9 years or by a 
rcdshift greater than 6 in the low density universe. 



in Newtonian collapse calculations as found by Ciotti et al. 
(2007). 

The cosmic time and redshift of the initial collapse of 
objects of different mass is shown in Fig. 4 as a function of 
mass. The initial conditions and final state of these objects is 
given in Table 1. The initial collapse is that point where the 
virial ratio IT /V reaches its first maximum corresponding to 
excessive kinetic energy primarily in tangential motion. We 
see that initial re-collapse for galaxy mass objects (10 1() — 
10 12 Mq) occurs almost coevally, when the cosmic age is 
between 50 and 100 million years, or at redshifts between 
10 and 20. By redshifts of 8 to 10, these galaxy scale masses 
would be in place as virialized objects. 

For determining the density and velocity distributions 
the effective number of shells is increased to 8800 by consid- 
ering the position and velocities of shells at different epochs. 
The density distribution resembles that of a Jaffe model 
(Jaffe 1983) having a power law of exponent near 2 in the in- 
ner region steepening to more than 3 in the outer regions (we 
recall that this is also similar to the density distribution in 
the MOND isothermal sphere: Milgrom 1984, Sanders 2000). 
The radial velocity dispersion gradually declines as a func- 
tion of radius, as in the high order polytropic spheres con- 
sidered by Sanders (2000). This is roughly consistent with 
decline observed in elliptical galaxies, at least, out to the 
effective radius. 

Because of the large collapse factor the final equilibrium 
form is extremely anisotropic with an effective anisotropy 
radius of 1 kpc- well within the effective radius. Such ob- 
jects would probably be highly unstable due to the ex- 
treme anisotropy, but this is suppressed here by the im- 
posed spherical symmetry. The non-linear development of 
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Figure 4. The rcdshift and cosmic age corresponding to initial 
re-collapse (defined as maximum virial ratio 2T/V) for masses 
ranging from 10 9 to 10 13 Mq. Galaxy scale objects all re-collapsc 
between redshifts of 10 and 20. 



this instability would lead to a triaxial figure and a larger 
anisotropy radius consistent with stability, as is the case in 
the three-dimensional collapse calculations of Nipoti, Lon- 
drillo & Ciotti (2007a). For this reason, the equilibrium ob- 
jects resulting from these spherically symmetric calculations 
are not appropriate in detail as models for elliptical galax- 
ies. Nonetheless, it is of interest that the global properties 
resemble those of the anisotropic polytropes applied previ- 
ously (Sanders 2000) as models for elliptical galaxies: specif- 
ically, a mild deviation from isothermality similar to that of 
n=12 to n=16 polytropes and a velocity distribution which 
is increasingly anisotropic in the outer regions. 

In Fig. 5 the surface density profiles are shown for the 
equilibrium figures resulting from four different calculations 
appropriate to different mass spheres. We see that these pro- 
files are well fit by empirically successful r 1//4 law typically 
over two orders of magnitude in surface density. The flat 
central surface density is an artifact of the unrealistically 
large initial tangential velocity on shells. Significantly, the 
projected central surface density is constant, and the effec- 
tive radius grows, roughly, as the square root of the mass. 
The effective radius is also compara ble to th e critical ra- 
dius for modified dynamics (r m = GM/ao)\ i.e., galaxy 
mass objects naturally re-collapse to a radius of about 10 
kpc without dissipation. 

It is of some interest that spherically symmetric expan- 
sion and re-collapse in MOND seems to be capable of pro- 
ducing final virialized objects in which the surface density 
distribution is a reasonable approximation to a r 1 ^ 4 law. 
With Newtonian dynamics, it is necessary to begin with a 
highly inhomogeneous, or clumpy, initial density distribu- 
tion in order to drive the redistribution in phase space suffi- 
cient to achieve such a universal surface density distribution; 
uniform spherical collapse does not work (van Albada 1982). 
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Table 1. Initial and final properties of spherical collapse calculations. (1): mass of 
sphere (10 11 M©; (2): initial cosmic time (current Hubble time); (3) initial radius 
in kpc; (4) initial velocity at outer radius (km/s); (5): redshift at re-collapse; (6) 
final los velocity dispersion (km/s); (7): final effective radius in kpc. 



Mass (10 11 M ) 
(1) 


tt (10- 3 Ho- 1 ) 
(2) 


r; (kpc) 
(3) 


Vmax (km/s) 

(4) 


^recoil 

(5) 


cto (km/s) 
(6) 


r e ff (kpc) 
(7) 


0.05 


1.66 


5.6 


100.2 


29.1 


72.1 


1.75 


0.10 


1.86 


7.0 


135.7 


20.7 


95.8 


3.75 


0.20 


2.10 


8.9 


163.9 


18.9 


108.0. 


5.25 


0.40 


2.37 


11.2 


198.0 


17.0 


139.4 


7.0 


0.60 


2.54 


12.8 


221.1 


16.6 


155.7 


8.5 


0.80 


2.68 


14.1 


239.1 


15.7 


172.7 


10.8 


1.0 
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Figure 5. The figures show the surface density distribution for 
four different MOND collapse simulations corresponding to the in- 
dicated masses. In all cases the simulation begins with the sphere 
in near Hubble expansion having an over-density of S = 0.01 
where the initial conditions are given by the numerical solution 
of eq. 6. The plots are of log surface brightness vs. r 25 and 
the dashed lines indicate the best r 25 law fit with the effective 
radius indicated. The equilibrium figures are reasonably approx- 
imated by the de Vaucouleurs law. The flattening in the central 
regions is artificial and due to the unrealistically large assumed 
tangential velocity on the spherical shells. The projected central 
surface brightness, however, is the same in each case. 



Of course, it is impossible to repeat the calculations de- 
scribed here with Newtonian dynamics because such spheres 
would be unbound; there would be no re-collapse. But be- 
ginning with cold spherically symmetric Newtonian collapse 
of objects of the same mass and size scale as described here, 



I find that the final range of phase space covered by the par- 
ticles is considerably more restricted than for the MOND 
calculations. The MOND re-collapse produces a number of 
shells which would be unbound in the context of Newtonian 
dynamics; this appears to make the critical difference for the 
final surface density distribution. 

Fig. 6 shows the results of the 13 simulations cov- 
ering a range of masses between 5 x 10 9 and 2 x 10 12 
Mq on the log effective radius-log central velocity dis- 
persion plane compared to the observations of j0rgensen 
et al.(j j0rgcnse n 1999 J0rgensen, Franx fc Kaergard 1 995a 



|j0rgensen, Franx fe Kaargard 1995b I. In several of these 
simulations a different value of the initial tangential veloc- 
ity was assumed (v t = 20-40 km/s). Again we see that these 
global characteristics of objects which condense out of the 
cosmological expansion via the MOND prescription are gen- 
erally similar to those of actual elliptical galaxies, although 
the models are somewhat more inflated than actual elliptical 
galaxies (larger effective radius for a given velocity disper- 
sion). This suggests that some global dissipation may be 
necessary to reproduce the observed distribution of real ob- 
jects on this plane. 

Fig. 7 is the mass- velocity dispersion relation for these 
13 equilibrium objects; also shown is that of the anisotropic 
n=12 to n=16 polytropes applied previously (Sanders 2000) 
as MOND models for elliptical galaxies which reproduce the 
observed fundamental plane. For the collapse models the 
velocity dispersion is that along a single line-of-sight toward 
the centre weighted by density. The observed line-of-sight 
velocity dispersion within a circular diaphragm, for example 
of 1.6 kpc diameter (|j0rgcnscn, Franx fc Kaargard 1995a[ ), 
would be about 25% lower. 

It is evident that the collapsed spheroids do exhibit a 
mass-velocity dispersion relation similar to that observed, 
although they are more homologous than actual ellipticals 
(less scatter on the Faber- Jackson relation) . This is because 
the spherical collapse model is highly idealised (i.e., no devi- 
ation from spherical symmetry or internal structure) . These 
objects would also lie on a fundamental plane similar to 
that described by the anisotropic polytropes (Sanders 2000, 
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Figure 6. The dark points show the distribution of the equi- 
librium collapse models on the log(cr )-log(r e ) plane. The cross 
points are observed elliptical galaxies from the samples of 
j0rgensen et al. The solid line is eq. 10 with p = 1. 



Figure 7. The dark points show the mass-velocity dispersion 
relation for the equilibrium collapse models; the fainter crosses are 
the same for the anisotropic MOND polytropes applied previously 
(Sanders 2000) as models of elliptical galaxies. The solid line is 
cq. 11 with g = l. 



eq. 15) i.e., Af/(lO n )M « 10~ V(km S - Y' 76 Mkpc)]°- 98 
(the scaling is lower due to the extreme and unrealistic 
anisotropy). But the important conclusion is that objects 
resembling the ellipticals do condense out of the Hubble flow 
within the framework of the MOND scenario for structure 
formation. Moreover, the final virialized objects attain the 
mean binding energy of actual galaxies without the necessity 
of dissipation. These are not deep MOND objects (such as 
dwarf spheroidals or other low-surface brightness systems) 
but are more similar to average massive ellipticals which are 
Newtonian within an effective radius. (Deep MOND isother- 
mal spheres with large constant density cores are possible 
configurations, but not by this formation scenario; Milgrom 
1984, and private communication.) 

In summary, objects formed from MOND dissipation- 
less collapse exhibit well-defined radius-velocity dispersion 
and mass- velocity dispersions relationships. These are of the 
form 



and 



r e = pa 2 /a 



qGMa , 



(10) 



(11) 



where p « q « 1 (shown by the solid lines in Figs. 6 and 7. 
These objects differ from pure isotropic, isothermal spheres 
where p = 4.36 and q = 0.0625 (Milgrom 1984, Sanders 
2000). Significantly, this means that such objects also have 
a characteristic density which is also related only to the 
velocity dispersion: 



Pe = 



2-Kqp z Ga 2 



(12) 



4 THE MOND CONDITION FOR 
FRAGMENTATION 

These dissipationless calculations still do not address the 
question of why galaxies of stars seem to be restricted to 
mass of less than 1O 12 M0. For this we have to consider cool- 
ing and fragmentation processes which lead to star forma- 
tion. 

In the traditional picture, this process is considered in 
terms of two competing timescales- the dynamical time vs. 
the cooling time. Presumably, the elements of the collapsing 
cloud collide, heat up, and attain a temperature appropri- 
ate to the virial velocity dispersion. The radiative cooling 
timescale of this hot plasma is 

t c = 3fcT[^m p A(r)n] _1 (13) 

where T is the temperature, p, is the mean molecular weight, 
m p the mass of the proton, n the density, and A(T) is the 
cooling rate per hydrogen nucleus in a plasma having spec- 
ified, presumably primordial, abundances. With Newtonian 
dynamics, the dynamical, or collapse, timescale is 



(14) 



where p is the density. 

Fig. 8 is a density-temperature plane for a hypo- 
thetical gas cloud and the solid curve shows the locus 
where t c ~ td- Above the solid curve, one would ex- 
pect cooling, fragmentation and star formation to play 
an important role. The points on the curve show three 
classes of pressure supported systems: globular clusters 
(Pryor & Meylen 1993, Trager et al. 1993), elliptical galax- 
ies Jj0rgensen 1999 |j0rgensen, Franx fc Kaargard 1995a[ 



|j0rgensen, Franx fc Kaargard 1995bJ and X-ray emitting 
clusters of galaxies (White, Jones & Forman 1997). The 
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Figure 8. The solid curve shows the defines the locus on the 
density-temperature plane where radiative cooling time equals 
the classical, Newtonian dynamical timescale. Gas clouds in the 
region above the curve can cool and fragment before the cloud 
collapses. Here one would expect to form self-gravitating objects 
consisting of stars. The parallel dotted lines are the locus of ho- 
mogeneous objects of indicated mass in virial equilibrium. The 
squares, crosses and triangles show, respectively, globular clus- 
ters, elliptical galaxies and X-ray emitting clusters of galaxies. 
The long dashed line is the locus of cooling time equal to Hubble 
time; cooling flow clusters should lie above this line. 



Figure 9. The solid curve shows the defines the locus on the 
density-temperature plane where radiative cooling time equals 
the MOND dynamical timescale. As above, gas clouds in the 
region above the curve can cool and fragment before the cloud 
collapses. Here one would expect to form self-gravitating objects 
consisting of stars. The vertical dotted lines arc the locus of ho- 
mogeneous objects of various mass in virial equilibrium in the 
context of MOND. As before, the squares, crosses and triangles 
show, respectively, globular clusters, elliptical galaxies and X-ray 
emitting clusters of galaxies. Note that the vertical scale differs 
from that of Fig. 8. 



point made by numerous authors (e.g., Silk 1977) is clear: 
Those objects consisting primarily of stars, ellipticals, are 
above the curve; whereas clusters of galaxies, consisting pri- 
marily of hot gas, are below. 

If we consider gas clouds uniformly mixed with a dark 
matter halo in virial equilibrium we can also define a relation 
between the mean density and the temperature for a given 
dark matter mass: 

P =¥h{— )V»itf- (15) 

where a is a number depending upon the density distribu- 
tion in the cloud and (ss 0.15) is the baryon to dark 
matter density ratio,. This density-temperature relation is 
also shown in Fig. 8 for objects of the indicated total dark 
mass. The point is that for objects with M > 10 12 Mq the 
virialized cloud will lie primarily in the region where cooling 
is slow compared to the collapse timescale. In this region 
we would expect that the fragmentation to the level of in- 
dividual stars does not occur before subsequent dynamical 
evolution of the object. 

With MOND we have seen that dissipationless objects 
which condense out of the Hubble flow exhibit a well defined 
radius-velocity dispersion, mass-velocity dispersion relations 
(eqs. 10 and 11) which imply a density- velocity relation (eq. 
12). However, the initial re-collapsing object, at maximum 
expansion, is certainly a gas cloud; as the cloud collapses 
subcomponents will collide and, most likely, generate a tem- 



perature with thermal velocity comparable to the random 
velocity of components. Therefore, a critical assumption 
here, as in the standard model, is that the gas virial temper- 
ature generated during re-collapse, fragmentation, and star 
formation is essentially given by the random velocity of the 
virialized object. For our hypothetical gas cloud, therefore, 
the radius-temperature and mass-temperature relations cor- 
responding to eqs. 10 and 11 are (for p = q = 1) 

re = 4 - 5 (iw) kpc (16) 

M / T \ 2 

i^= ai4 y 

where I have taken the mean atomic weight of 0.62 for 
an ionised gas with primordial abundances. The dynami- 
cal timescale for the cloud would therefore be t d — pa/ao 
or 

* d = 3.8 x 10 7 (-^) 2 years (18) 

Expressing the characteristic density of a collapsed object 
(eq. 12) in terms of the hydrogen density for a fully ionised 
plasma with primordial abundances: 

nff =°- 9 (iwr lcrrr3 (w) 

The condition t c = td where now td is denned by eq. 
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18 is shown by the solid curve in Fig. 9. The three classes 
of self-gravitating objects are are indicated here as is Fig. 8. 
The temperature of the gas clouds with indicated masses, 
as given by eq. 17, is now shown by the vertical dotted 
lines (there is no dependence on scale, and hence density, 
in the MOND virial relation). And finally, the characteris- 
tic MOND density (eq. 19) for a near-isothermal cloud in 
equilibrium is shown by the dashed line. 

We see that the solid curve again separates those objects 
built primarily from stars (elliptical galaxies) from objects 
consisting primarily of gas (clusters). We also see that the 
MOND density-temperature relation intersects this curve 
for an object with a mass of about 10 12 Mq. This means 
that objects which are more massive will have a density and 
temperature such that cooling and fragmentation occurs on 
a timescale longer than subsequent evolution of the system 
as a whole; this provides a natural upper limit for objects 
built of stars. It is also evident that actual virialized systems 
do cluster about the line of the characteristic MOND den- 
sity. This has been pointed out previously in other contexts 
(Sanders & McGaugh 2002). 

If cosmological neutrinos have a mass of 2 eV they will 
comprise the dark component of clusters required even in 
the context of MOND (Sanders 2003). In this case the neu- 
trino fluid will become a dominant component of clusters at 
temperatures higher than about 4 x 10 7 K, and the curve 
of td = t c would steepen at higher temperatures (Sanders 
2007). This does not affect the argument on fragmenta- 
tion. If the dark component of clusters is baryonic (Milgrom 
2007), then the form of the curve remains the same. 

At high temperature, the solid curve in Fig. 9 (t c — t d ) 
appears to approach a constant value of the density. This is 
because the cooling timescale for a high temperature plasma 
is set by free-free emission and is given by 

tc = (20) 
/i 2 n e m e 2 rt 2 c 2 a 

where K is a numerical constant, m e is the mass of the 
electron, r t is the Thompson radius (e 2 /m e c 2 ) and a is the 
fine structure constant (e 2 /Tic). That is to say, the cool- 
ing timescale is proportional to VT, as is the dynamical 
timescale (eq. 18). Therefore, the temperature drops out of 
the condition t c = id in this high temperature regime; there 
is a characteristic density where this condition is met of 
n w 0.1 cm -3 . Equating this to the MOND density for a 
virialized system (eq. 19) then defines a characteristic tem- 
perature or (from eq. 17) a mass. This mass would be the 
upper limit to an object where fragmentation to stars has 
proceeded- the characteristic galaxy mass. This turns out 
to be 

M c «a 6 a 9 3 (^) 3 ^«lO 12 M (21) 
V m e / c 3 

where a g — Grn p 2 /he is the gravitational coupling constant. 
Therefore, with MOND, as in the standard paradigm, the 
galaxy mass is defined in terms of fundamental constants 
although, in our case, ao enters into the expression. Given 
that, coincidentally, Tiolq/c 3 w 0.1a 9 ra e , eq. 15 becomes 

M c ^0.1a 6 a 9 - 2 (^) 2 m p (22) 



which is similar (but not identical) to the expression derived 
by Silk (1977). 



5 CONCLUSIONS 

It has been appreciated for many years (e.g. Binney 1977) 
that, by the mechanism of gravitational instability in an 
expanding Universe, it is impossible to form the luminous 
parts of galaxies without dissipation. This remains true in 
the context of the current cosmological paradigm, ACDM. 
Given the amplitude of fluctuations on galaxy scale at de- 
coupling, bound objects of 10 12 Mq would form at a red- 
shift of about 5 with a characteristic scale of to 20 to 50 
kpc (defined as the radius at which the power law density 
distribution breaks from -1 to -3) and a density between 
10~ 24 and 10~ 25 gem -3 - independently of whether or not 
they form hierarchically or by monolithic collapse. (Navarro, 
Frenk & White 1997). But this would be the presumed char- 
acteristic density of dark matter. If the baryonic matter were 
distributed in the same way (no dissipation) then its den- 
sity would be 10 times lower. The baryonic matter must 
dissipate, radiate its kinetic energy and collapse by a fur- 
ther factor of 10 to produce the observed baryonic density 
in the inner luminous regions of galaxies. Moreover, this dis- 
sipation must be global; fragmentation to the level of stars 
before collapse would not produce galaxies of the observed 
surface brightness. 

With MOND, given the simple ansatz that the modified 
gravity only applies to fluctuations, this is not true. Then, 
in a pure baryonic Universe, galaxy scale fluctuations sepa- 
rate out of the Hubble flow early (z > 35) and re-collapse 
and virialize by a redshift of 10. The mean surface density 
and binding energy of these collapsed objects is comparable 
to that of observed galaxies. The standard surface bright- 
ness distribution for ellipticals, the r 1//4 law, is reproduced 
as well as the Faber-Jackson and Fundamental Plane cor- 
relations. This is all accomplished without dissipation. Of 
course, galaxies are self-gravitating ensembles of stars. So, to 
form galaxies in this dissipationless limit, cooling and frag- 
mentation to the level of stellar mass objects must occur on 
a timescale short compared to the collapse time. Therefore, 
as in the standard picture, radiative cooling must set the 
mass scale of galaxies. 

The more robust result here, independent of the ansatz 
on structure formation, concerns the properties of proto- 
galactic clouds. With MOND, unlike the Newtonian case, an 
acceleration scale enters the structure equation. This means 
that a near isothermal object with a given velocity disper- 
sion, or temperature, will have definite mass: if that velocity 
dispersion is 100-200 km/s the mass will be on the order of 
10 11 Mq Moreover, there also is a definite size scale asso- 
ciated with an isothermal object of a given mass- and that 
is the radius of the transition from Newtonian dynamics to 
modified dynamics. The isothermal object will inflate to this 
radius and thereafter truncate. This implies the existence of 
a characteristic density which can be identified with a partic- 
ular mass; for example, a mass of 10 11 Mq will have a density 
of about one particle per cubic centimetre. For lower masses, 
or temperatures below 10 7 K, the virialized MOND clouds 
lie in the domain where cooling is more rapid than collapse; 
here we expect the cloud to fragment and form a galaxy of 



© 0000 RAS, MNRAS 000, 000-000 



10 R.H. Sanders 



stars essentially by dissipationless collapse. For higher tem- 
peratures this is no longer the case. The cloud may well 
maintain its identity as a gaseous object until merging with 
a similar cloud as part of a larger collapsing structure. It is 
of interest that observed pressure supported, near isother- 
mal objects do, at present, lie near the characteristic MOND 
density-temperature relation. 

In summary, with MOND galaxy scale masses are likely 
to re-collapse and virialize early (z > 10). Spherically sym- 
metric N-body calculations indicate that the objects which 
condense out of the Hubble flow, resemble actual elliptical 
galaxies. This suggests that elliptical galaxies may form by 
monolithic collapse without global dissipation. The condi- 
tion that cooling takes place rapidly compared to collapse, 
as in the standard scenario, places an upper limit of about 
1O 12 M0 on those bound objects which can consists primarily 
of stars. 

I am grateful to Francoise Combes, Stacy McGaugh, 
and Moti Milgrom for very helpful comments on this paper. 
I also thank the referee, Pasquale Londrillo, for a number of 
very usful comments and criticisms which greatly improved 
the content and presentation of this paper. 
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